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ABSTRACT 

We present a detailed description of the first direct measurement of the spatial correlation function 
of X-ray selected AGN. This result is based on an X-ray flux-limited sample of 219 AGN discovered 
in the contiguous 80.7 deg^ region of the ROSAT North Ecliptic Pole (NEP) Survey. Clustering is 
detected at the 4ct level at comoving scales in the interval r = 5 — 60 Mpc. Fitting the data 
with a power law of slope 7 = 1.8, we find a correlation length of r-Q = 7.4^i['g Mpc {Qm = 0.3, 
ft\ = 0.7). The median redshift of the AGN contributing to the signal is — 0.22. This clustering 
amplitude implies that X-ray selected AGN are spatially distributed in a manner similar to that of 
optically selected AGN. Furthermore, the ROSAT NEP determination establishes the local behavior of 
AGN clustering, a regime which is poorly sampled in general. Combined with high-redshift measures 
from optical studies, the ROSAT NEP results argue that the AGN correlation strength essentially 
does not evolve with redshift, at least out to z ^ 2.2. In the local Universe, X-ray selected AGN 
appear to be unbiased relative to galaxies and the inferred X-ray bias parameter is near unity, 6x 1- 
Hence X-ray selected AGN closely trace the underlying mass distribution. The ROSAT NEP AGN 
catalog, presented here, features complete optical identifications and spectroscopic redshifts. The 
median redshift, X-ray flux, and X-ray luminosity are z = 0.41, fx = 1.1 x 10"^'^ ergs cm~^ s~^, and 
Lx = 9.2 X 10-*^ 

ergs s ^ (0.5-2.0 keV), respectively. Unobscured, type 1 AGN are the dominant 
constituents (90%) of this soft X-ray selected sample of AGN. 

Subject headings: quasars: general — cosmology: observations — large-scale structure of universe — 
X-rays: general 



I 1. INTRODUCTION 

Active galactic nuclei (AGN) are presumably the very 
luminous manifestations of accretion onto supermassive 
black holes. Shining brightly across the electromagnetic 
spectrum and easily detectable to very high redshift, 
AGN are accessible tracers of galaxy formation and evo- 
lution as well as large-scale structure. The spatial distri- 
bution of AGN reflects the distribution of matter fluctua- 
tions modulated by the complex, non-linear astrophysics 
of black hole formation. Thus measurements of AGN 
clustering and its evolution provide important tests for 
models of AGN formation in an adopted cosmological 
reference frame ijHartwick fc Schad3ll990() . 

The vast majority of work on this topic has focused 
on optical surveys for quasi-stellar objects (QSOs). 
Following the first attempt to measure AGN cluster- 
ing by Osmer (1981), t he first significant detection 
came from the work of iShaveil (^984). In the en- 
suing years sample si zes have increased as has the 
result i ng precis ion (Sh anks et al.lll987Hlovino fc Shavep 
19881 lAndreani fc Cristianil [199 2t 'Mo fc Fang' '1993; 
Shanks Rovlfj 119941 \cZnm k Shanks 1991 
La Franca et al.l Il998j) . There is broad agreement 
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in the conclusions of these studies. Clustering is de- 
tected at the ~3-4(T level. The spatial distribution is 
characterized using the conventional approach of the 
two-point correlation function which tracks the excess 
probability over random to find two objects separated 
by a given distance. This function is found to have a 
power-law shape with a slope of 7 ~ 1.8. The clustering 
scale length, where the correlation function is unity, 
is ro ^ 6 h^^ Mpc measured at a mean redshift of 
z ~ 1.4. This amplitude is comparable to that found for 
luminous, local galaxies, and larger than that of dwarf 
or low surface brightness gala xies (e.g., the resu lts from 
the Sloan Digital Sky Survey: Zchavi et al. 2002 and the 
2dF Galaxy Redshift Survey: Hawkins ct al. 2003). This 
is consistent with the idea that moderately powerful 
AGN are found in moderately luminous galaxies, and so 
share their clustering properties. 

One question left unresolved by the aforementioned pa- 
pers is the nature of the clustering evolution , the redshift- 
depen dent behavior of the scale length. iCroom et alJ 
(l2ml have used the 2dF QSO Redshift Survey (2QZ) 
to definitively redress this issue, at least in the high red- 
shift interval, and have found essentially no change in 
the clustering amplitude between z ~ 0.5 and z ~ 2.2 
(see also iC room et al. 2003). The recent measure by 
iGrazian et al.i (^20041 of the local AGN population lends 
further support to this largely flat behavior. 

At X-ray wavelengths there are far fewer results on the 
clustering of AGN chiefly due to a historic lack of suffi- 
ciently large X-ray samples, particularly those with com- 
prehensive optical follow-up. The i?0S'j4T all-sky survey 
(1990-1991), the first with an imaging X-ray detector, 
and the subsequent program of pointed observations en- 
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abled the first significant steps forward. Despite the cur- 
rently limited development, this field is anticipated to 
experience rapid growth with the advent of much deeper 
surveys from Chandra, XMM-Newton and future X-ray 
survey missions. 

The use of X-rays in the selection and characterization 
of AGN is a natural choice. X-r ay emission appe ars to 
be a univer sal feature of AGN ijElvis et all (1971) . and 
essentially all optically selec ted AGN are X-ray luminous 
l|Avni fc Tananbaumlll98 6ll. In fact, X-ray emission is 
likely the least biased selecti on technique, par ticularly 
at hard energies above 2 keV l|Mushotzkvll2f)0^ . 

There have been several determinations of the angu- 
lar correlation of X-ray AGN (Vikhlinin & Forman 1995; 
Akvlas et ani200(i iGiaccon i et alji200H iBasilakos et all 
20041) . However, translating these results to constraints 
on the three-dimensional clustering strength requires a 
number of assumptions that results in only weak con- 
straints on tq. Actually measuring the spatial correla- 
tion function is more involved since optical identifications 
and redsh ift measures are required. Early a ttempts by 
iBovle fc M o (1993) and lCarrera et all (|1998() using data 
from Einstein and ROSAT, respectively, did not yield 
significant clustering signals. 

The first direct measure of the spatial correlation func- 
tion of X-ray selected AGN was achieved in the context 
of the ROSAT North Echptic P ole (NEP) Sur vey. This 
result was initially reported by iMullisI l|200H) . In the 
current paper we provide a detailed description of this 
clustering analysis. Note that only very recently have 
additional measures of this kind beco me availabl e as re - 
ported for the Chandra deep fields bv lGilli et al.' f20G4'). 

We describe the ROSAT NEP Survey in §|51 focusing 
on the basic properties of the X-ray flux-limited AGN 
sample and the survey selection function. In §|31 the 
formalism of the clustering analysis is laid out, includ- 
ing a description of the spatial correlation estimator and 
the Monte Carlo simulations of the survey volume. We 
present the spatial correlation function of X-ray selected 
AGN in §0] These results are compared to similar stud- 
ies at both X-ray and optical wavelengths with emphasis 
on the evolution of clustering strength with redshift in 
§|S1 We close with a summary of the key results in §|S1 

Unless otherwise stated, we use a cosmological model 
with the parameters. Ho = 70 /170 km Mpc~^, 
= 0.3, and fl\ — 0.7, and refer to this as ACDM. 
However, for the correlation length we use Hq = 100 h~^ 
km s~^ Mpc~^ for ready comparison to previous results. 
Unabsorbed X-ray fluxes and restframe X-ray luminosi- 
ties are quoted in the 0.5-2.0 keV energy band, unless 
otherwise indicated. Measurement errors are given at the 
68% confidence interval (la). 

2. THE ROSAT NEP AGN SAMPLE 
In the ROSAT All-Sky Survey (RASS), an 80.7 
deg2 region around the NEP (aaooo = IS'^OO™, 
^2000 = -l-66°33') constitutes one of the deepest ob- 
servations of the X-ray sky e ver achieved with such a 
large , contiguous solid angle ()Mulh£ll2001l : IHenrv et al.l 
I200I IVoges ei~al]l200l . Here 445 unique X-ray sources 
are detected with fluxes measured at greater than 4(t 
significance in the 0.5-2.0 keV energy band. The median 
and maximum exposure times are approximately 5 and 
38 ks, respectively. We have identified the physical 



nature of 443 (99.6%) of the ROSAT NEP X-ray sources 
through a c omprehensive pr ogram of imaging and 
spectroscopy (|Gioia et al.ll2003|) . AGN are the dominant 
constituents comprising nearly half (49.4%) of identified 
sources in the ROSAT NEP Survey catalog (Figure Pi . 

The ROSAT NEP AGN are particularly well suited for 
a clustering analysis. The AGN are drawn from a con- 
tiguous, wide-angle region in the sky sampled to a rela- 
tively deep X-ray flux. The survey selection function is 
well determined and is only a function of X-ray flux. We 
have identified all but two of the 445 X-ray sources in the 
survey region. Thus the AGN sample is essentially com- 
plete and requires no complicated assumptions to correct 
for incompleteness. Furthermore, we have spectroscopi- 
cally measured redshifts for the entire sample. 

2.1. Basic Properties 

We present the fuU sample of 219 ROSAT NEP 
AGN in Tabled This includes several notable re visions 
relative to previous versions of the catalog (|Mullij 
2001; Gioiactal. 2003). First, the AGN fluxes and 
luminosities previously reported were over-estimated by 
approximately 20% on average due to an error in the 
conversion of X-ray count rate to flux.^ Second, the 
sample has grown by one due to the re-classification 
of an X-ray source (RX J1824. 7+6509; see footnote in 
Table nj. And finally, we have adopted the presently- 
favored "concordance" cosmology in computing X-ray 
luminosities. The revised and updated catalog with 
corrected properties is presented in this paper and 
should be the reference point in any future work with 
the ROSAT NEP AGN sample. Below we explain the 
columns of Table 

Columns (1) and (2): the object name and internal 
identification number. Sources are listed in order of in- 
creasing right ascension. 

Columns (3) and (4): the right ascension and declina- 
tion of the X-ray centroid, respectively. 

Columns (5) and (6): the right ascension and declina- 
tion of the optical counterpart, respectively. 

Column( 7 ): Col umn density of Galactic hydroge n from 
|eI vis et al.l ((1994j) with supplements from .Stark et alJ 
II1992D . 

Columns (8) and (9): Count rate and count rate error 
in the 0.1-2.4 keV energy band measured within a cir- 
cular aperture of 5' radius. The quoted error is the la 
uncertainty based on a maximum-likelihood analysis. 

Column (10): Unabsorbed flux in the 0.5-2.0 keV 
energy band derived from the count rate assuming the 
source has a power- law spectrum (photon index F = 2) , 
with absorption fixed at the Galactic value for the source 
position. This total flux includes a correction factor of 
1.0369 which accounts for the small fraction of point 
source flux falling outside the photometry aperture. 

Column (11): Rest-frame, K-corrected luminosity in 
the 0.5-2.0 keV energy band. Note the K-correction for 
a power-law spectrum is (1 + z)^^'^ and is thus unity 
for r = 2 (cosmology: Hq — 70 /170 km s^^ Mpc~^, 
= 0.3, and = 0.7). 

^ A power-law photon index of r=l was mistakenly used instead 
of r=2 due to a typographical error in the analysis software. X-ray 
fluxes and luminosities for galaxy clusters and stars in the refer- 
enced works are not affected by this problem. 
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Fig. 1. — Distribution on the sky of the 219 AGN in the ROSAT NEP Survey (Aitoff projection). To give some perception of the 
three-dimensional arrangement of the sample, the symbol size reflects the projected size of a constant physical radius of 500 kpc at 
the redshift of the AGN. As the X-ray emitting region of AGN is sub-parsec, this scale is greatly exaggerated for visualization purposes. 
An animated "fly-through" of the ROSAT NEP Survey volume is available at http : //www. if a.hawaii . 6du/~ mullis/nepSd.html. 



Column (12): Spectroscopic redshift (typical uncer- 
tainty < 0.001). 

Column (13): AGN classification based on the equiv- 
alent width of emission lines and the broadness of 
permitted lines (1 — type 1,2 = type 2). 

We have secured spectroscopic redshifts for the en- 
tire ROSAT NEP AGN sample; the fuU range is z = 
0.026 — 3.89 and the median is z = 0.41. X-ray lumi- 
nosities range over 10^^ — 10^^ ergs s^^ with a median 



of Lx = 9.2 X 10*3 g^gg g^i gj^^g ROSAT has very 
limited sensitivity to hard X-rays above 2 keV, the bulk 
of the NEP AGN (90.4%, 198 objects) are type I AGN 
(QSOs and Seyfert 1 galaxies) based on the equivalent 
width and the broadness of their permitted emission lines 
{Wx > 5 A, FWHM > 2000 km s^^). The remaining 21 
NEP AGN (9.6%) are categorized as type II (Seyfert 2 
and star- forming galaxies). 

2.2. Sky Coverage and logN(>S)-logS Distribution 
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Fig. 2. — Map of the X-ray flux limits for the detection of AGN in the ROSAT NEP Survey (Aitoff projection). The logarithmicaUy- 
spaced contours indicate the minimum total flux that an object must have to meet the 4(t detection requirement of the survey. Contours 
are labeled in units of lO"^** ergs cm~^ s~^ in the 0.5-2.0 kcV energy band. The enhanced sensitivity toward the center of the survey 
region is a result of the ROSAT scan pattern that converges at the NEP. 



The distribution on the sky of the 219 ROSAT NEP 
AGN is shown in Figure ^ The average AGN source 
density is 2.7 deg^'^ and increases toward the NEP at the 
center of the survey region as a resuh of the RASS scan 
pattern. During the survey phase of the mission, ROSAT 
scanned the sky in great circles which overlapped at the 
ecliptic poles. The resulting peak in integrated expo- 
sure time at the NEP is reflected in the sensitivity map 
shown in Figure El This map of limiting X-ray flux is 
derived from RASS exposu re and background maps (see 
Chapter 3 of iMulhsl I2OOII and J. P. Henry et al. 2004, 
in preparation). Here the contours, in units of 10^^** 
ergs cm~^ s~^, indicate the minimum flux that an AGN 
must have to produce at least a 4(t detection and thus 
meet the selection criterion of the ROSAT NEP survey. 
The sensitivity map is integrated to obtain the final se- 
lection function shown in Figure 01 The effective sky cov- 
erage is 81 deg^ at bright fluxes, starts to decrease below 
1.3 X 10"^'^ ergs cm~^ s~^, and remains significant down 
to ^ 2.3 X 10^^^ ergs cm^^ s^^ where the coverage is 1 
deg^. 



In subsequent modeling of the AGN population 
we will make use of the integral number counts or 
logiV(> S") -log5 distribution of the ROSAT NEP 
AGN sample. We calculate the number of objects per 
square degree observed above a flux S by summing up 
the contribution of each source, weighted by the area in 
which the source could have been detected, via the equa- 
tion, 

where Si is the flux of the «th source and ^{Si) is the 
sky coverage at flux 5*^ f rom Figure (HI A maximum likeli- 
hood fit of a power law ({Murdoch. Crawford, fc .Tauncevl 
119731) to the differential source counts over the flux range 
(2-1000) X 10^^** ergs cm^^ s~^gives the relation 

/ q \ -l-30±0.08 

iV(>5) = (1.5±0.1)(^^^j deg-^ (2) 

The normalization is set by the total number of AGN 
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observed above the minimum survey flux. The mea- 
sured log A^(> S) — logiS* and the power-law fit are plot- 
ted in Figure^ These results are consistent with pre- 
vious determin a tions in this flux regime. For example, 
iHasinger et"al] (|199,1D found a slope of 1 72 ±0.27 at 
fluxes > 2.7 X 10"^'' ergs cm~^ s~^, and iMason et al.l 
slope of 1.6 ±0.3 at fluxes above 
ergs cm~^ s~"^. 



I|2000() report 
3 X 10-" 



3. CLUSTERING ANALYSIS 

One of the simplest and most popular techniques for 
characterizing the clusterin g of objects is the two-point 
spatial correlation function (iPeeblesllQSOj) . This statistic 
is defined in terms of the joint probability (dP) of simul- 
taneously finding an object in a volume dVi and another 
object in a volume dV2 separated by a distance r, 



dP = 



ar)] dVidV2, 



(3) 



where n is the average number density of objects and (r ) 
is the two-point spatial correlation function. In a uniform 
random Poisson point process, objects are distributed 
in a completely random pattern and the probabilities of 
finding objects in dVi and dV2 are independent such that 
^(r) = 0. If objects are more clustered than average 
^(r) > 0, whereas, if objects are more dispersed than 
average ^(r) < 0. Hence, the two-point correlation is 
the excess probability over random to find two objects 
separated by distance r. 

The observed correlation function for a variety of ex- 
tragalactic objects (e.g., galaxies, AGN, and galaxy clus- 
ters) is well fit by a power law of the form 



(4) 



where the slope is typically 7 « 1.8 while the correlation 
length ro depends on the type of object. 

3.1. Correlation Estimator 

The correlation function is extracted from spatial data 
via a pairwise analysis. In essence, the number of ob- 
ject pairs of a given spatial separation in the data are 
compared to the number of object pairs of the same sep- 
aration from a random catalog. The random catalog is 
created by homogeneously populating the survey volume 
in a manner that is consistent with the selection func- 
tion and boundary conditions of the survey. Hence, an 
enhancement of the number of data-data pairs relative 
to the corresponding random-random pairs is indicative 
of structure in the data. 

Though the correlation function can be com- 
puted using the natu ral form, ^ri(r) ~ (DD/RR) — 1, 
iLandv fc Szalavl l)1993t) have demonstrated that the vari- 
ance is minimized using the estimator. 



CM 



N,{N, - 1) DD{r) (N, - 1) DR{r) 



N^Nd-l) RR{r) 



Nd RR{r) 



1. (5) 



Here DD{r) is the number of data-data pairs in the NEP 
AGN sample with separations of r ± Ar/2 in redshift 
space. Similarly, RR{r) is the number of random-random 
pairs and DR{r) is the number of data-random cross- 
pairs, each with separations of r ± Ar/2. The number 
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Fig. 3.— Sky coverage for AGN in the ROSAT NEP Survey. 
This is the total area of the sky surveyed as a function of total flux. 
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Fig. 4. — Integral number counts per square degree of AGN 
in the ROSAT NEP Survey. The dashed line is a power-law fit 
to the data derived from a maximum likelihood treatment of the 
differential number counts. The shaded region is the la error region 
for this fit. 



of objects in the data sample is A^d, and the number 
of objects in the random sample is A^r- A number of 
different cor relation estimators are com monly used (e.g., 
IDavis fc Peebles 1983; Hewctt 1982; Hamilton'199 3), but 
a deta iled comparison by Kerscher. Szapudi. fc Szalavl 
()2000() shows that the Landy & Szalay estimator is pre- 
ferred in terms of minimal bias and variance, and gives 
the most reliable results on all spatial scales. 

We calculate co-moving separations r between ob- 
jects in redshift space using the standard r elations (e.g., 
IWeinbergllT97l IPeeblesI [T993t iHoggl IT991 . The trans- 
verse co-moving distance between the observer and an 
object at redshift z is 



Dm{z) 



5 / dz'[QMil + z'f + 



(6) 
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where ^Im is the matter density parameter, and Q/^ is the 
cosmological constant. The curvature of space is charac- 
terized by f2fc = 1 — — n\. The function S is defined 
according to: 



S{x) 




ilfc > open universe 
rife = flat universe 
rife < closed universe 



Note that the transverse co- moving distance (Dm) is re- 
lated to the luminosity distance and the angular di- 
ameter distance (Da) via the equations _Dl = (1 + z)Dm 
and Da = Dm/(1 + z). 

Consider two objects separated on the sky by an angle 
9 with transverse co-moving radial distances of I?mi and 
Dm2- The CO- moving separation of the second object as 
measured from the first object is 




1)^2 - 2 d Dmi Dm2 cos ( 



(7) 



where 




(8) 



fe.g.. lOsmeilll98lHMatarrese et al.lll997|) . In a flat uni- 
verse rife = 0, d = 1 and Equation[7|reduces to the cosine 
rule for Euclidean space. Note that Equation is not 
symmetric for the positions of the two objects. However, 
this is only important in the case in which d is not close 
to unity. 

3.2. Construction of Random Samples 

Generating random samples for the survey volume is 
critical for estimating the correlation function. The 
objective is to construct simulated datasets with con- 
stituents that are consistent with the physical distribu- 
tion of the source population and are selected in the same 
manner as the real data. In the following Monte Carlo 
procedure we create an X-ray flux-limited sample of sim- 
ulated AGN possessing the same flux and redshift dis- 
tributions as the actual ROSAT NEP AGN. Of course 
the important difference between the simulated and real 
AGN is that the parent population of the former is as- 
sumed to be randomly distributed on the sky and in red- 
shift space. 

To properly model the non-uniform sensitivity pattern 
of the ROSAT NEP survey region (Figure 0), the simu- 
lated AGN must have the same flux distribution as the 
real AGN. We have demonstrated in S I2.2I that the ob- 
served AGN flux distribution follows a power law of the 
form A^(> S) = KS^°' where a « 1.3 (FigurelH). Thus 
the differential probability distribution of fluxes scales 
like S'-l^+i). In practice the required set of random 
numbers with a power-law distribution is obtained from 
a set of random uniform numbers using a transformation 
method (e.g., Bcvington & Robinson 1992). If p is a ran- 
dom number uniformly distributed between and 1, then 
a random sampling in X-ray flux S above a limiting flux 
S'lini is distributed like 



Fig. 5.— Redshift distribution of the ROSAT NEP AGN. The 
black line is the distribution smoothed with a Gaussian of width 
<Jz = 0.1. The median redshift of the sample is 2 = 0.41. To 
preserve the readability of this plot, the most distant AGN at 
z = 3.889 is not shown. 



S = S'lini (1 -P) 



-1/q 



(9) 



The limiting flux in the ROSAT NEP Survey for AGN 
is S'lini = 2.0 X 10""'^'' erg s""'^ cm~^. Small variations to 
the value of a do not change the outcome of the clustering 
analysis. 

Redshifts for the simulated AGN are drawn from a 
probability density function based on the observed dis- 
tribution of the ROSAT NEP AGN (Figure EJ). Here the 
redshift distribution has been smoothed with a Gaussian 
of width (Tz =0.1. This smoothing scale is selected such 
that the resulting distribution has the overall shape of the 
data but is not strongly affected by discrete large-scale 
structure or Poisson fluctuations. Our results are insen- 
sitive to reasonable variations in this smoothing length. 

We randomly place simulated AGN within the ROSAT 
NEP Survey boundaries (17''15" < azooo < 18'H5™, 
62° < ^2000 < 71°). For each object we assign an X-ray 
flux according to Equation IHl and a redshift drawn from 
the smoothed distribution. Then we test to see if its 
flux is above the local flux limit for its particular posi- 
tion in the survey region as shown in Figure 12 Objects 
above the threshold are retained, while those below are 
rejected. We proceed to randomly populate the survey 
volume in this manner until we have constructed a ran- 
dom catalog of 100,000 simulated AGN; i.e., ~ 450 times 
the size of the ROSAT NEP AGN sample. 

4. RESULTS 

We present the spatial correlation function of the 
ROSAT NEP AGN in Figure El Here the data and sim- 
ulated objects have been binned up in terms of their pair 
separations (e.g., DD{r), RR{r), DR{r)) and evaluated 
using the adopted correlation estimator (Equation EJ. 
The uncertainties associated with these data points are 
typically estimated using Poisson statistics of the form. 
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(10) 
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Fig. 6. — Spatial correlation function of the ROSAT NEP AGN. The error bars on the data points are la Poisson errors. The dashed 
line is the maximum-likelihood fit (ro = 7.4^j'g Mpc, 7 = 1.8 fixed) to the range from 5 to 60 Mpc {i^M = 0-3 and Ha = 0.7). 
The shaded region demonstrates the ilcr errors on the correlation strength. 



where DD{r) is the number of data-data pairs in the in- 
terval. However, the Poisson distribution is not well ap- 
proximated by a Gaussian in situations where the num- 
ber of counts in each bin is small {N < 20) and thus the 
^/N errors in the denominator of Equation llOl can under- 
estimate the 68 % confid e nce le vel. To avoid this we use 
the formulas of iGehreT^ l)1986j) to estimate the Poisson 
confidence intervals for one-sided 84% upper and lower 
bounds which correspond to ±la in Gaussian statistics. 

A positive clustering signal is readily apparent in the 
spatial correlation function of ROSAT NEP AGN and is 
significant at the 4ct level based on the excess AGN pairs 
relative to random at r < 60 Mpc. Thus it is desir- 
able to estimate the correlation strength and slope using 
the canonical power-law fit. This could be done using 
the coarsely binned data in Figure |3 and minimizing the 
statistic. However, this approach is subject to uncer- 
tainties because the determination can be sensitive to the 
size and distribution of the selected bins. The maximum- 
likelihood method is an alternative means for determin- 
ing a power-law fit to the correlation function which has 
the advantage of making maximal us e of the data and 
is free from arbitrary binning (e.g., |(,^r9f1^ ^1^ ^-lIll^^Tt. 
Popowski et al. | [ 19981: Borgani. Plionis. fc Kolokotroniif 
~ l200nHCollins et'a l. 2000). The 
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co-moving separation r is parsed into very small inter- 
vals such that there is either or 1 object pair in any 
given interval. In this limit Poisson probabilities are ap- 
propriate. The probability P of observing Vi object pairs 



where fii pairs are expected is given by 



(11) 



In the sparse sampling limit, the probabilities associated 
with the bins are independent of each other such that a 
likelihood function C can be defined in terms of the joint 
probabilities. 



n 



This leads to the useful expression 



In/: 



+ i^i In/ii - In j/J) . 



(12) 



(13) 



Here the summation is over all the intervals within 
the range of co-moving separations for which the power- 
law fit is to be determined. The number of observed 
object pairs i^i are the data-data pairs measured from 
the ROSAT NEP AGN, whereas the expected numbers 
of pairs fii are calculated using the Landy & Szalay es- 
timator in Equation El Notice that ^i, which is DD{ri) 
in Equation[3 is a function of the power-law parameters 
ro and 7 and the pairwise data from the observed and 
random samples, specifically DR{ri) and RR{ri). 

The best-fitting values of the power-law fit to the cor- 
relation function are determined by minimizing the ex- 
pression 
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5 = -2 InC (14) 

with the confidence levels defined to be 

AS* = 5'(ro,7) - 5(ro,bcst,7bcst)- (15) 

Since S is distributed like x^, the la, 2a, and 3a (68.3%, 
95.4%, and 99.7%) confidence intervals for a two p aram- 
eter f i t are AS = 2.30, 6.17, and 11.8, respectively l|Avnil 
ITOTfit ICasi3IT976l |T979). If we fix the slope, the confi- 
dence intervals for the one-parameter fit are AS* = 1.0, 
4.0, and 9.0. 

We show the results of the maximum-likelihood analy- 
sis of the ROSAT NEP AGN clustering data in FigureEl 
The power- law fit over the range from 5 to 60 Mpc 
indicates best-fit values of rg = 7.5^42 Mpc and 
7 = 1.851q gQ where the errors are la for two interest- 
ing parameters. If we set the slope at 7 = 1.8 (the value 
typically found for normal galaxies and optically selected 
AGN), the correlation length is tq = 7.4ti g h'^ Mpc. 
We overplot this best-fit power-law on the pairwise data 
in Figure El The lower limit on the fitting range is set 
by the smallest AGN- AGN separation in the NEP data; 
the upper limit matches the general break in the power- 
law shape of the observed function. Our results do not 
change significantly for reasonable variations in the fit- 
ting range. For example, varying the upper limit between 
40 Mpc and 80 h^^ Mpc results in the best-fit ro for 
7 = 1.8 ranging between 7.3 Mpc and 7.6 Mpc 
(i.e., Arg < 0.2a). Note that the effective redshift of 
the NEP clustering signal is = 0.22 given by the me- 
dian redshift of the AGN in pairs with separations in 
the range 5-60 ft.^^ Mpc. This characteristic redshift is 
smaller than the sample median (z=0.41) since the den- 
sity of detected objects, and thus the number of close 
pairs, decreases with redshift in this flux- limited survey. 

Repeating this analysis assuming an Einstein-de Sit- 
ter cosmology (JIm = 1, = 0), we find best-fit 
values of tq = 7.2tl-^h-^ Mpc and 7 = 1.83j:}:[;J], and 
ro = 7.1 ± 1.6 h^^ Mpc with 7 = 1.8. This modest de- 
crease in the correlation length is expected since the ef- 
fect of a significant flA is to increase object-object sep- 
arations r and thus the resulting correlation length rp. 
These results are in excellent agreement with the prelim- 
inary analysis of the NEP data first reported in Mullis 
ll200l . 

There is a substantial supercluster of galaxies 
at z = 0.087 in the ROSAT NEP Survey region 
ijMullis et alJl200'l|l . This is manifested as 5a spikes in 
the redshift distributions of galaxy clusters and IRAS 
galaxies, but only a 2a fiuctuation in AGN density. 
There are 12 AGN from our sample in the redshift regime 
of the NEP supercluster (0.07 < z < 0.1). These objects 
comprise 5% of the total sample and 12% of the AGN 
at z < 0.4. We have examined the potential impact of 
the superstructure on the clustering signal by restricting 
the analysis to the 200 AGN at z > 0.1. This result is 
entirely consistent with the results from the analysis of 
the fuU sample (Aro ~ 0.6cr, A7 ~ 0.2a). Thus the NEP 
supercluster has no significant effect on our results. 

A final concern lies in the potential consequences of 
source confusion as a result of the angular resolution 
of the ROSAT PSPC detector. The FWHM of the 
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Fig. 7. — Probability contours for the power-law normaliza- 
tion (ro) and slope (7) from a maximum-likelihood analysis of 
the ROSAT NEP AGN (Q^/ = 0.3 and Qa = 0.7). The soUd 
contours with shading are the la, 2a, and 3a (68.3%, 95.4%, 
and 99.7%) confidence intervals for a one parameter fit. The 
dashed contours are the corresponding intervals for a two parame- 
ter fit. The filled circle marks the best-fit for two free parameters, 
ro = 7.5^^ 2 ^"""^ Mpc and 7 = LSS^g gp; the open circle marks 
the best- fit for one free parameter with the slope fixed at 7 = 1.8, 
ro =7.4^j'g?t~^ Mpc. The dotted horizontal line indicates the 
canonical power-law slope for the correlation function (7 1.8). 



RASS PSF is ^35", and thus we expect sources sepa- 
rated by less than ^--^30" to be unresolved. However, 
the confusion-limited regime is never approached due to 
the relative sparsity of AGN at the sensitivity limits of 
the ROSAT NEP Survey. For example, the smallest dis- 
tance between two NEP AGN is 3.8'. Furthermore, at the 
typical redshift of AGN contributing to the correlation 
detection (z^ = 0.22), 30" corresponds to a co-moving 
separation of 56 h^^ kpc which is much smaller than the 
best-fit correlation length. Thus confusion effects should 
not bias the correlation analysis. 

5. DISCUSSION 

In this section we discuss the significance of the 
ROSAT NEP results for the clustering of AGN, make 
comparisons to similar measures, and discuss the impli- 
cations of our results for the evolution of AGN clustering. 

5.1. First Direct Measure of the Spatial Correlation 
Function of X-ray Selected AGN 

Our results from the analysis of ROSAT NEP data 
represent the first direct and significant measure of the 
spatial correlation function of X-ray selected AGN. This 
outcome is due to the advantageous combination of rel- 
atively deep sensitivity over a wide, contiguous survey 
region which results in sufficient numbers of AGN- AGN 
pairs at small scales to yield a significant signal. For ex- 
ample, the maximal signal-to-noise ratio in the pairwise 
data is observed at co-moving separations of less than 
30 h-^ Mpc. Here we find 123 AGN- AGN pairs where 
only 78 pairs are expected if the data had a uniform spa- 
tial distribution. Thus AGN clustering is detected at 
greater than 4a level. 

The correlation length derived from the ROSAT NEP 
AGN sample is consistent with the value ro ^ 6 Mpc 
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assoc iated with normal galaxies (e.g., |H^^m^^^ 
1 200.'^ and optically selected AGN (e.g.. iCroom et a" 
1200 1*). This suggests that X-ray luminous AGN are spa- 
tially clustered in a manner similar to that of the afore- 
mentioned objects. In fact it lends additional evidence to 
the idea that X-ray selected and optically selected AGN 
are drawn from the same population. Our observations 
imply that AGN are not biased relative to galaxies based 
on the similarity of clustering strengths. Note that the 
factor of two density enhancement observed in both AGN 
and galaxies within the NEP supercluster provi des fur- 
ther support to this view (see Figure 1 in iMuUis et al.1 
120011) . As observation and theory i ndicate the bias of 
local galaxies is close to unity fe.g.. iVerde et alJl2002t 
iWeinberg et al.l 120041 and references therein), the in- 
ferred bias parameter for X-ray selected AGN is also 
near unity, 6x ~ 1- Previous measurements of the X-ray 
bias, some of which inc luded all extragalact ic sources, 
have varied widely (see IBarcons et al.ll2'001l and refer- 
ences therein) ; however our NEP measurement is in good 
agreement w ith recent results base d on the hard X-ray 
background ijBoughn fc Crittendenl l2"004). 

It is pertinent to review the two efforts to measure 
the spatial correlation of X-ray selected AGN prior to 
that of the NEP. The fact that there are such few stud- 
ies reflects the difficulty of constructing sufffciently large 
and appropriately distributed samples to make this mea- 
sure. The first clus tering analysis of this type was per- 
formed bv lBovle fc Mo (1993) using 183 AGN (z < 0.2) 
taken from the Einstein Extended Medium Sensitivity 
Survey (EMSS; Stocke et al. 1991). Though this sample 
is similar in size to the NEP, it is spread over a nearly 
ten times larger solid angle (770 deg^). Moreover, their 
survey region is not contiguous, but rather the combi- 
nation of hundreds of pointings distributed across the 
sky. [BovIc & Mo (1993) found the significance of the 
clustering signal was only O.Scr and concluded that AGN 
clustering at small scales is very weak, if present at all. 

The second work to pursue the spa tial correlation func- 
tion is that of ICarrera et alJ (jl998|) who examined two 
separate samples of X-ray selected AGN. The first con- 
sisted of 107 objects (0 < z < 3.5, z = 1.43) fro m the 
ROSAT Deep Survev fPRS: IShanks et^HmS) . The 
second was composed of 128 AGN (0 < z < 3.5, z = 
0.84) from th e ROSAT Internati onal X-ray Optical Sur- 
vey (RIXOS: iMason et alJl2000|) . Like the EMSS, the 
solid angles surveyed by the DRS and RIXOS (1.4 deg^ 
and 20 deg^, respectively) are th e summation o f many 
individual pointed observations. ICarrera et al.l (|1998f) 
found no indication of clustering in the DRS sample, but 
detected clustering at r < 40-80 h^^ Mpc in the RIXOS 
sample at about 1.7(7 significance. Given the limitations 
of the data, the authors did not present an explicit cor- 
relation function. However, by combining the DRS and 
RIXOS AGN, they derived a correlation length in the 
range of 1.5 < ro < 5.5 Mpc dependent upon the 
assumed model of clustering evolution. These results for 
an Einstein-de Sitter cosmological model would increase 
by approximately 30% in the concordance model. 

Aside from our ROSAT NEP results, the only other 
significant measure of this kind to date is from the recent 
analysis of the Chandra deep fields. IGilli et al] (j2004f) 
have measured the projected correlation function (i/;(rp) 
where Vp is the co-moving separation perpendicular to 



the line of sight) using the 2Msec Chandra Deep Field 
North (CDFN) and IMsec Chandra Deep Field South 
(CDFS). This work is based on 160 AGN from the CDFN 
(mean redshift, z - 0.96) and 97 AGN from the CDFS 
(z ~ 0.84). Though these samples are smaller than that 
of the NEP, they are concentrate d in much s maller solid 
angles, ^ 0.1 deg^ each. Thus [gTiH et al.1 are able to 
extract high signal-noise results by probing the correla- 
tion function at very small scales (~ 0.2 — 10 Mpc) 
where the signal is comparatively stronger. Converting 
the projected correlation to the three-dimension scale 
length assuming a power-law shape, they find for the 
CDFN ro = 5.5 ± 0.6 h-^ Mpc, 7 = 1.50 ± 0.12; for the 
CDFS ro = 10.3 ± 1.7 Mpc, 7 = 1.33 ± 0.14. If the 
slope is fixed at their preferred value of 7 = 1.4, they find 
for the CDFN tq = 5Ato.t Mpc and for the CDFS 
ro = 10.4 ±0.8 h^^ Mpc. There is a large variance be- 
tween the two fields. However the strong enhancement 
in the correlation strength of the CDFS is attributed to 
two large redshift spikes at z ~ 0.7. 

It is enlightening to compare the ROSAT NEP and 
CDFN/S spatial correlation functions. The clustering 
strength of the NEP, ro ~ 7.4 Mpc, lies intermedi- 
ate in the ro = 5 — 10 Mpc range of the Chandra deep 
fields. The ROSAT PJG^ are a soft X-ray (0.1-2.4keV) 
selected sample compared to the Chandra samples as- 
sembled from detections at both soft to hard energies 
(0.5-10 keV). Consequently the i?OS'ylT NEP sample is 
dominated by type 1 AGN, whereas those of CDFN/S 
are composed of nearly equal mixes of type 1 and type 2 
AGN. However . .Gilli et al.. (,,2004|) have measured the cor- 
relation strength separately for hard versus soft detected 
AGN as well as type 1 versus type 2 AGN and find 
no significant differences. In relative terms, the NEP 
AGN are luminous and nearby while the CDFN/S AGN 
are faint and distant. The NEP characterizes the local 
(z ~ 0.2) clustering of AGN versus the distant (z ~ 0.9) 
measures of the CDFN/S, a point we will elaborate on 
in The median luminosity of the ROSAT NEP 

AGN is ix = 9.2 X lO''^ h^^ ergs s~^ whereas the mean 
of the CDFN/S is Lx = 6.0 x 10''^ /i'^ ergs s'^ (both in 
the 0.5-2.0 keV energy band). 

5.2. Comparison to the Angular Clustering of X-ray 
Sources 

Thus far we have restricted our discussion to the direct 
measure of the spatial correlation of AGN using X-ray 
flux- limited samples. There are very few investigations 
of this kind due to the observational challenges of secur- 
ing X-ray data of sufficient depth and/or breadth, and 
obtaining optical follow-up observations for object iden- 
tifications and spectroscopic redshifts. If we allow for ad- 
ditional assumptions and uncertainties, the angular cor- 
relation function of X-ray sources can be incorporated. 

The angular clustering o f X-ray sources has been 
detected several ti mes (e.g., Vikhlinin & Forman 19951 
Akvlas et alJl200(it IG iacconi et al. 2001: Basilakos et all 
20041) . The bulk of this signal is presumably generated by 
AGN since they are the dominant class of X-ray emitters 
at high Galactic latitudes. However, in the absence of op- 
tical follow-up, the exact AGN fraction is uncertain. The 
completely identified ROSAT NEP Survey provides a ref- 
erence point at soft energies (0.5-2.0 keV) and at mod- 
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Fig. 8. — The scale length of AGN clustering as a function of lookback time (lower axis) and redshift (upper axis). X-ray results are 
shaded symbols; optic al results are solid symbols. Symbol key and referenc es for the data are: AERQS (filled square; RJrazian ct al. 200 j), 
RASS (shaded circle; Akylas ct al.' 2000), NEP (shaded star; MuUis et aLll2001l and this paper), CDFN/S Cshaded triangles: iGilh et all 
120041) . and 2QZ (filled circles; Groom ct al. 2003). All of these are from direct measures of t he spatial correlatio n function except for the 
RASS. Over-plotted are theoretical predictions based on different biasing models described bv lGroom et al.l 120011 and references therein): 
linear theory (solid line), long-lived AGN (dotted line), merger model for different values of the minimum halo mass (dashed line, lO^^ Mq, 
10^'' Mq, 10^^ Mq bottom to top). The dot-dashed line is the best-fitting empirical bias model for the 2QZ data. The uncertainty intervals 
of the measurements are based on one-parameter fits. The NEP and RASS use a fixed slope of 7 = 1.8; AERQS and 2QZ assume 7 = 1.56. 
CDFN/S use their preferred value of 7 = 1.4. All of the measurements and models are for a AGDM cosmological model. 



erately faint fluxes (a few times 10 ergs cm ^ s 



where ~50% of the X-ray sources are AGN. A second 
important caveat concerning these angular studies is 
that both the redshift distribution of the sources and 
the redshift dependence of the clustering evolution must 
be assumed to extract the three-dimensional correlation 

le ngth rp. 

iVikhlinin fc FormanI l)1995j) reported the first detec- 
tion of the angular clustering of discrete X-ray sources 
using a set of deep ROSAT observations. They con- 
cluded that the correlation strength is consistent with 
that of optically selected AGN assuming such objects 
constitute a l arge fraction (>50%) of their X-ray sources. 
lAkvlas et al.l l|2000ll measured the angular clustering of 
X-ray sources from the RASS Bright Source Catalog, tak- 
ing precautions to exclude stellar and extended sources 
as much as possible. Assuming a typical redshift of 
z ~ 0.1—0.2, comoving clustering evolution, and a power- 
law slope of 7 w 1.8, they derived a correlation length 
of ro « 6.4 ± 1.7 Mpc; thus similar to the results at 
high redshift from optical surveys. We have converted 
their EdS result to AGDM assuming a median redshift 

of z = 0.15.^ 

Recently iBasilakos et al.l l)2004fl analyzed the cluster- 
ing of XMM-Newton hard X-ray sources (2-8 keV) in 



^) a 2 deg^ region. Exploring a variety of potential lu- 



minosity functions (and thus redshift distributions) and 
clustering characteristics, they found a rather high cor- 
relation length span ning the range rn ~ 9 — 19 
Mpc for 7 = 1.8. IBasilakos et alJ argue that hard- 
band sources are more strongly correlated than soft-band 
sources. This claim is supported by the counts-in-cell 
analysis from deep Chandra data by Yang et al. (2003), 
but is at odds with the results of iGilli et al.t (^2004) for 
the CDFN/S. 

5.3. Evolution of AGN Clustering 

To round out this discussion, we incorporate the latest 
results on the clustering of optically selected AGN, and 
use the ensemble of X-ray and optical data to constrain 
the evolution of AGN clustering. 

As outlined in §n there is a significant body of work in 
the literature characterizing the clustering of optically se- 
lected AGN. T he current sta t e of th e art is encapsulated 
in the work of iCroom et"an 1)200 IL HOOS) with the 2QZ 
survey. Analyzing a sample of over 20,000 AGN, they 
find the two-point spatial correlation has a power-law 
shape over the range 1 — 60 Mpc. Their best-fit pa- 
rameters are tq = 5.76lo;27 h'^ Mpc and 7 = 1.64t^;[J^ 
(AGDM) measured at a mean redshift of redshift z ~ 1.5. 
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Moreover, the 2QZ sample is sufRciently large to map 
out the varia t ion in clustering strength at high redshifts. 
iCroom et al.l l|2003j) find no evolution in the clustering 
amplitude from z ~ 0.5 to z 2.2. 

Ironic in the era of high-redshift studies, the clustering 
properti es of local AGN has onlv recentlv been measured 
directly (IMulhsllmnt IGrazian et al.ll200l . Our ROSAT 
NEP analysis has an effective redshift of ~ 0.2. 
IGrazian et alJ lj2004|) have determined a similar measure 
for optically selected AGN using data from the Asiago- 
ESO/RASS QSO Survey (AERQS). They find an ampli- 
tude of ro = 8.6 ± 2.0 hr^ Mpc (ACDM) at an effective 
redshift of = 0.06. The AERQS and NEP results are 
in excellent agreement. 

A composite view of the clustering evolution of AGN 
is shown in Figure |H1 Here we plot the scale length of the 
spatial correlation as a function of redshift and lookback 
time; X-ray results are shaded symbols, optical results 
are solid symbols. Note that all of these are direct mea- 
sures of the spatial correlation function except for the 
RASS data point (Akvlas et al. 2000) which is a trans- 
formation of the angular signal. Any conclusions to be 
drawn from these data must be taken with due regard for 
several caveats that will be discussed shortly. However, 
for now we take the data points at face value. 

To the zeroth order no evolution is observed in the 
clustering strength of AGN out to redshift z ~ 2.2. 
A straight line at ro = 6 ft.^^ Mpc is a good fit to 
the data with only two points deviating by greater 
than 3(7. The CDFS is an obvious outlier, but, as 
noted before this high value is likely due to cosmic 
variance. Looking with more detail, a weighted mean 
of the local measures (NEP, RASS, and AERQS) indi- 
cates ro = 7.3 ± 1.0 Mpc; whereas the high-redshift 
results (2QZ and CDFN/S) have a weighted mean of 
ro = 5.80 ± 0.21 h'^ Mpc. Thus there is an -25% in- 
crease in clustering at low redshift [z ^ 0.2) relative to 
high redshift (z « 1 — 2). However, this very mild redshift 
evolution is significant at only the XAa level. 

The conclusion of zero to weak negative evolution is 
generally applicable to either the X-ray or optically se- 
lected AGN on their own, though with much less con- 
fidence for the former. Our ROSAT NEP results and 
the Chandra deep fields provide complementary redshift 
coverage. However, the trend based on strictly these X- 
ray AGN is still open ended due to the high variance of 
CDFN/S results. Forthcoming analysis at high redshift 
from Chandra and XMM-Newton will hopefully remove 
this uncertainty. We anticipate future findings will favor 
correlation lengths of ro « 6 Mpc on average. 

Our interpretation of the clustering behavior shown in 
Figure IHl must be tempered by several potentially impor- 
tant issues, including: 1) X-ray versus optical selection of 
the AGN, 2) comparison of different AGN types, 3) com- 
parison of different AGN luminosities, and 4) comparison 
of amplitudes derived for different power-law slopes. The 
first two concerns are strongly coupled. Both soft X-ray 
selection and optical selection preferentially detect the 
unobscured, type 1 AGN. Conversely, hard X-ray selec- 
tion also recovers the obscured, type 2 AGN. Assuming 
the unified paradigm for AGN is true (i.e., type 1 ver- 
sus 2 properties are mainly due to viewing geometry), 
there is no obvious reason to expect the two classes to 
cluster differently. This anticipation is supported by the 



tests of lGilh et all pOO^ in the CDFN/S, though the 
st rong angular clusterin g of hard X-rav s ources claimed 
bv lYang et all l)2003() and lBasilakos et alJ l|2004D could be 
weak counter examples. Note that the veracity of AGN 
unification is debated, and arguments can be made for 
the i ndependent evolution of t he two AGN populations 
(e.g., iFranceschini et af] 120021 and references therein). 
Except for the deviant CDFS result. X-ray and optical 
AGN cluster similarly at similar redshifts based on the 
data in Figure |H1 

The third issue of comparing dissimilar luminosi- 
ties is relevant if this property is correlated with the 
mass of the dark matter halo in which AGN re- 
side. Many popular clustering models make this as- 
sumption though its validity is not yet observation- 
ally confirmed (e.g.. Groom et al. 2002). First consider 
the median luminosities of the X-ray selected samples. 
The ROSAT NEP AGN, with a median luminosity of 
Lx = 9.2 X 10"*^ h^Q ergs s~^, are on average about fif- 
teen times more luminous than those from the CDFN/S 
(Lx = 6.0 X lO'*^ ergs s^^). We estimate the me- 
dian luminosity of the RASS AGN oflAkYlas et al. ( 2000) 
to be a few times lO'*'^, thus intermediate to the NEP 
and CDFN/S. It is interesting to note that the median 
X-ray luminosity of the ROSAT NEP AGN is rather sim- 
ilar to those that we estimate for the 2QZ and AERQS. 
Taking the mean absolute magnitudes of these two sam- 
ples (2QZ: Mb = -25.1 1 . ICroom et a].ll20fil AERQS: 
Mb = -22.99. ICrazian et al.ll2004l both Hq = 70 km 
s^^ Mpc^^, ACDM) and using an aver age spectral en- 
ergy distribution fe.g.. lElvis et alJll994D or the observed 
X-ray/optical relation (e.g.. lVignali et ani2003(l . we find 
mean luminosities of Lx ~ 4 x lO'^^ h^^ ergs for the 
2QZ and Lx « 1 x lO^'' /i^ ergs s'^ for the AERQS. 

The last complication to consider in regards to the data 
in Figure IHl lies in the fact that different authors use dif- 
ferent values for the power-law slope of the correlation 
function in deriving the clustering amplitude. The NEP 
and RASS analyses, for example, are based on a fixed 
slope of 7 = 1.8. The 2QZ and AERQS adopt a slope 
of 7 = 1.56. The CDFN/S report single-parameter fits 
for their preferred value of 7 = 1.4. Thus problems may 
arise in comparing the clustering strengths since the scale 
length and the slope are correlated to some degree. How- 
ever, the magnitude of the effect appears to be manage- 
able. For instance, in the case of the NEP these parame- 
ters are not strongly linked (see Figure [TJ . The NEP cor- 
relation length decreases by ~14% to ro ~ 6.5 h^^ Mpc 
if we lower the slope to 7 = 1.56. Similarly, Gilli et ^ 
(2004) report their best-fit correlation lengths for the 
CDFN/S each increase by 15% if they assume a slope 
of 7 ~ 1.8. Notice such adjustments to homogenize the 
surveys tend to erase the marginal evolution trend dis- 
cussed above, and further strengthen the case for no red- 
shift evolution. 

The above caveats notwithstanding, we briefly out- 
line possible theoretical interpretations of the clustering 
strength behavior as a function of redshift. The man- 
ner in which mass clusters via gravitational instability is 
generally well understood. Thus the spatial correlation 
function of matter fluctua tions, £ma,ss. can be computed 
with fittle ambiguity fe.g.. iSmith et alJl2003() . However, 
connecting this mass clustering to AGN clustering is non- 
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trivial since we do not know the distribution of AGN 
relative to the distribution of mass. This bias between 
the luminous and dark matter densities is parameterized 
by b{z) and reflects the unknown physical mechanisms of 
AGN formation. 

The evolution of the AGN correlation function is often 
formulated in the following terms, 

C(r, z) = b^z) D\z) U.ssir, z = 0). (16) 

Here D{z) is the linear growth of density perturbations 
which is {l + z)~^ for an EdS cosmology. The present-day 
spatial correlation function of mass is ■^mass('', ^ = 0). 

We introduce three models for consideration. The sim- 
plest scenario of clustering evolution assumes the bias 
does not change with redshift. The linear theory model 
scales like D'^{z) and is represented by the solid line in 
Figure IHl The second model, the long-lived scenario, pro- 
poses that AGN are formed at an arbitrarily high redshift 
and then move thereafter according to the gravi tational 
potential of the density fluctuations l)Frvlll99^ . Here 
the bias parameter increases with redshift which coun- 
teracts the linear growth factor and produces negative 
evolution which is less steep than linear theory (dotted 
line in Figure |HJ). The third biasing scenario, known as 
the merger model, uses the Press-Schechter l)1974l) for- 
malism to track the evolution of the dark matte r halo 
mass function (Matarrese ct al. 1997; Moscardin i et al.l 
11998ft . The resulting bias parameter evolves more steeply 
than the growth of perturbations enabling a model clus- 
tering strength that increases with redshift. The over- 
all strength of clustering reflects the mass of the dark 
matter halos hosting the AGN. The merger model pre- 
dictions for several minimum halo masses are plotted in 
Figure IHl (dashed lines). More ext ensive discuss i ons o f 
thes e biasing models are laid out in iCroom et al.l l)200lD 
and lGrazian et alJ l|2n04D . 

The linear model is clearly too steep and inconsistent 
with the flat to mildly negative evolution delineated by 
the clustering observations. The less steep long-lived 
model agrees with the local measures but can not repli- 
cate the flat behavior of the 2QZ data at high redshift. In 
terms of overall shape, the closest agreement is given by 
the merger m odel. The empirical biasing description of 
ICroom et al.l (HOQIJ is based on the merger model. Their 
best fit with a minimum halo mass of ~10^^ M© is 
shown as the dot-dashed line in Figure |S| The ROSAT 
NEP clustering strength is consistent with this model as 
are most of the other data, except for the CDFS. 

6. CONCLUSIONS 

We have made the first direct measurement of the spa- 
tial correlation of X-ray selected AGN. These results are 
based on the analysis of 219 AGN detected at soft ener- 
gies (0.1-2.4 keV) in the ROSAT NEP Survey. The AGN 
catalog, presented here, includes updated X-ray and op- 
tical properties and features complete optical identifica- 
tions and spectroscopic redshifts. The clustering signal 
is significant at the 4ct level corresponding to a cluster- 



ing length of To = 7.4j]i[ g h^^ Mpc assuming a power-law 
shape with a slope 7 — 1.8. The median redshift of the 
AGN contributing to this signal is = 0.22. 

Our results indicate that X-ray and optically selected 
AGN share similar clustering properties as both display 
clustering scale lengths of rg ~ 6 Mpc. This is not 
too surprising given that soft X-ray and optical selection 
are preferentially sensitive to unobscured type 1 AGN. 
Assuming that we are probing the same population, the 
ROSAT NEP measure provides an effective zero-point 
in the study of clustering evolution. This low redshift 
determination is quite valuable since optical surveys are 
largely ineffective in this regime. The no-evolution trend 
delineated by the 2QZ at high redshifts appears to extend 
to low redshifts based on the NEP analysis. We find a 
modest increase in clustering strength (^25%) at 2; < 0.2 
relative to the high redshift results at z ^ 1.5; however, 
this is only 1.4cr significant. 

Since the amplitudes of galaxy and AGN clustering 
are similar, we argue that AGN are not biased relative 
to galaxies. It likely that AGN randomly sample the 
galaxy distribution and do not preferentially probe the 
high peaks of the matter density field. Given that the 
bias of local galaxies is near unity, the inferred X-ray 
bias is near unity. Hence it seems that X-ray selected 
AGN closely trace the underlying mass distribution. 

The study of X-ray AGN clustering is still in the 
early stages and much work lies ahead. The NEP 
has provided the first firm data point, and the recent 
results from the CDFN/S provide small uncertainty but 
large variance measures at high redshift. Forthcoming 
results from on-going and planned deep surveys with 
Chandra and XMM-Newton will play a vital role by 
hopefully averaging out cosmic variance and securing 
the clustering amplitude at high redshift. Future X-ray 
survey missions like the Dark Universe Observatory will 
provide precision measures of X-ray AGN clustering as a 
function of high redshift with the ROSAT NEP results 
serving as an important anchor point at low redshift. 
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0.24 


0.3533 


1 




RX J1745. 7+6748 


2610 


■IT Ar An f 

17 45 42.6 


1 T AO '\ v 

+67 48 15 


IT Ar AC\ A 

17 45 42.4 


1 /^T AQ '\ A 

+67 48 14 


4.28 


0.0058 


0.0013 


5.99 


0.37 


0.4143 


2 




KX J 1745.9+6451 


2650 


17 45 55.2 


+64 51 18 


17 45 55.5 


+64 51 25 


3.49 


0.0211 


0.0024 


19.70 


0.18 


0.1790 


1 




RX J1746. 0+6727 


2670 


17 46 03.0 


+67 27 09 


17 46 01.8 


+67 27 09 


4.70 


0.0120 


0.0015 


12.94 


0.17 


0.2146 


1 




IV j'V d 1 i ^4:0 . 1 -j- U i O 1 


97nn 






J- / Uo.o 




4 21 


1464 




1 zlQ QS 


n OR 

U.VJU 




\ 




RXJ1746.2+6227 


2710 


17 46 14.6 


+62 27 01 


17 46 13.9 


+62 26 54 


3.25 


0.0319 


0.0034 


28.73 


413. ,35 


3.8890 


1 




RXJ1746.3+6320 


2750 


17 46 21.6 


+63 20 06 


17 46 21.8 


+63 20 10 


3.02 


0.0165 


0.0026 


14.31 


0.67 


0.,3697 


1 




RXJ1747.0+6836 


2800 


17 47 00.3 


+68 36 26 


17 46 59.9 


+68 36 34 


4.42 


0.2131 


0.0060 


223.35 


0.21 


0.06,30 


1 




RXJ1747.1+6813 


2810 


17 47 10.6 


+68 13 19 


17 47 12.7 


+68 13 26 


4.54 


0.0057 


0.0013 


6.05 


27.06 


2.3920 


1 




RXJ1747.2+6532 


2820 


17 47 14.4 


+65 32 30 


17 47 13.9 


+65 32 35 


3.66 


0.0108 


0.0016 


10.33 


15.11 


1.5166 


1 




RXJ1747.3+6702 


2840 


17 47 22.2 


+67 02 06 


17 47 21.5 


+67 02 01 


4.51 


0.0066 


0.0011 


6.99 


1.75 


0.7421 


1 




RXJ1747.4+6626 


2850 


17 47 26.8 


+66 26 27 


17 47 27.0 


+66 26 24 


3.78 


0.0138 


0.0014 


13.42 


0.07 


0.1391 


1 




RXJ1747.4+6924 


2860 


17 47 27.0 


+69 24 55 


17 47 27.9 


+69 25 09 


3.49 


0.0099 


0.0020 


9.25 


1.02 


0.5292 


2 
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Object NEP ax <5x aopt 5opt nn Count Rate Count Rate /x ^x ^ AGN 

ID (J2000) (J2000) (J2000) (J2000) [lO^O cm-^] [s^l] Error [s-l] [lO-^" cgs] [10"''' cgs] Type 



y'-) 


\'^) 








(0) 


(6) 




{' ) 


^^^^ 




{lU) 




(LZ) 






RX J1747.9+6538 


2890 


17 47 58.0 


+65 38 


35 


17 47 57.9 


+65 


38 


28 


3.85 


0.0305 


0.0022 


29.91 


1.04 


0.3248 


1 




RXJ1748.2+7016 


2900 


17 48 17.4 


+70 16 


14 


17 48 19.6 


+70 


16 


09 


3.93 


0.0446 


0.0036 


44.20 


0.43 


0.1858 


1 




RX J1748.3+6403 


2910 


17 48 22.7 


+64 03 


27 


17 48 23.1 


+64 


03 


38 


3.27 


0.0155 


0.0023 


14.00 


7.06 


0.9859 


1 




RX J1748.6+6842 


2940 


17 48 38.8 


+68 42 


11 


17 48 38.3 


+68 


42 


17 


4.25 


0.0236 


0.0023 


24.28 


0.02 


0.0537 


1 




RXJ1749.3+6411 


3000 


17 49 20.4 


+64 11 


08 


17 49 19.5 


+64 


11 


19 


3.27 


0.0109 


0.0019 


9.85 


4.94 


0.9836 


1 




RXJ1749.7+6422 


3001 


17 49 42.4 


+64 22 


46 


17 49 44.1 


+64 


22 


58 


3.24 


0.0078 


0.0017 


7.02 


1.83 


0.7540 


1 




RXJ1750.2+6814 


3050 


17 50 14.3 


+68 14 


33 


17 50 16.1 


+68 


14 


37 


4.50 


0.0069 


0.0013 


7.29 


0.12 


0.2310 


1 




RXJ1750.2+6415 


3060 


17 50 15.5 


+64 15 


15 


17 50 15.1 


+64 


14 


56 


3.11 


0.0115 


0.0019 


10.13 


0.19 


0.2504 


2 




RXJ1751.0+6710 


3100 


17 51 02.4 


+67 10 


09 


17 51 01.2 


+67 


10 


14 


4.33 


0.0038 


0.0008 


3.95 


0.56 


0.5870 


1 




RXJ1751. 1+6753 


3120 


17 51 09.5 


+67 53 


07 


17 51 08.9 


+67 


53 


08 


4.45 


0.0062 


0.0011 


6.52 


0.25 


0.3406 


1 




RXJ1751. 6+6540 


3160 


17 51 39.7 


+65 40 


40 


17 51 36.9 


+65 


40 


30 


4.10 


0.0122 


0.0015 


12.33 


4.02 


0.8259 


1 




RXJ1751. 9+6551 


3190 


17 51 57.6 


+65 51 


20 


17 51 56.7 


+65 


51 


17 


4.13 


0.0121 


0.0012 


12.28 


0.65 


0.3901 


1 


O 


RXJ1752.2+6624 


3210 


17 52 12.6 


+66 24 


56 


17 52 11.7 


+66 


24 


54 


3.88 


0.0049 


0.0008 


4.82 


0.27 


0.4002 


1 


RXJ1752.9+6440 


3230 


17 52 57.4 


+64 40 


58 


17 52 56.9 


+64 


40 


56 


3.22 


0.0125 


0.0016 


11.21 


0.04 


0.1230 


1 


a 


RXJ1753.1+6746 


40 


17 53 09.7 


+67 46 


44 


17 53 09.6 


+67 


46 


32 


4.77 


0.0070 


0.0011 


7.60 


5.37 


1.1297 


1 




RXJ1753.5+6811 


3231 


17 53 30.9 


+68 11 


47 


17 53 32.4 


+68 


12 


01 


5.01 


0.0056 


0.0014 


6.21 


0.43 


0.4366 


2 


lERING 


RXJ1753.6+6542 


3260 


17 53 41.6 


+65 42 


42 


17 53 42.1 


+65 


42 


40 


4.34 


0.0076 


0.0009 


7.90 


0.04 


0.1400 


1 


RXJ1753.9+7016 


3280 


17 53 55.5 


+70 16 


47 


17 53 56.6 


+70 


16 


42 


4.25 


0.0095 


0.0018 


9.77 


0.01 


0.0620 


1 


RXJ1754.0+6613 


60 


17 54 05.4 


+66 13 


54 


17 54 04.8 


+66 


13 


50 


4.12 


0.0152 


0.0010 


15.40 


0.90 


0.4067 


1 


RXJ1754.7+6819 


3340 


17 54 42.3 


+68 19 


08 


17 54 42.0 


+68 


19 


06 


5.23 


0.0066 


0.0011 


7.46 


0.27 


0.3292 


1 


RXJ1754. 7+6208 


3350 


17 54 43.2 


+62 08 


21 


17 54 42.3 


+62 


08 


30 


3.41 


0.0147 


0.0024 


13.58 


0.45 


0.3190 


1 


O 


RXJ1754.8+6706 


100 


17 54 49.3 


+67 06 


00 


17 54 49.7 


+67 


05 


56 


4.53 


0.0023 


0.0005 


2.44 


1.03 


0.9190 


1 




RX J1755.0+6446 


3360 


17 55 00.0 


+64 46 


32 


17 55 00.8 


+64 


46 


32 


3.49 


0.0085 


0.0013 


7.94 


1.65 


0.6870 


1 




RX J1755.0+6235 


3361 


17 55 03.6 


+62 35 


30 


17 55 03.2 


+62 


35 


41 


3.33 


0.0083 


0.0020 


7.57 


13.86 


1.6607 


1 


pi 


RX J1755.0+6519 


3370 


17 55 05.8 


+65 19 


50 


17 55 05.6 


+65 


19 


55 


3.78 


0.0948 


0.0030 


92.17 


0.14 


0.0785 


1 




RXJ1755. 1+6719 


110 


17 55 09.0 


+67 19 


50 


17 55 08.3 


+67 


19 


54 


4.69 


0.0049 


0.0008 


5.28 


0.08 


0.2225 


2 




RX ,11755.1+6852 


3380 


17 55 11.9 


+68 52 


30 


17 55 10.7 


+68 


52 


34 


4.56 


0.0062 


0.0012 


6.59 


7.43 


1.3645 


1 




RXJ1755.6+6209 


3420 


17 55 40.3 


+62 09 


41 


17 55 40.3 


+62 


09 


39 


3.41 


0.0240 


0.0027 


22.17 


0.04 


0.0846 


2 


H 


RXJ1755.6+7009 


3410 


17 55 40.5 


+70 09 


52 


17 55 41.3 


+70 


09 


51 


4.25 


0.0077 


0.0016 


7.92 


0.53 


0.4295 


1 


t-i 
H 


RX J1755.7+6249 


3430 


17 55 46.2 


+62 49 


27 


17 55 45.9 


+62 


49 


29 


3.33 


0.0402 


0.0034 


36.67 


0.61 


0.2360 


1 


O 


RXJ1755.9+6540 


140 


17 55 56.9 


+65 40 


54 


17 55 56.8 


+65 


40 


52 


4.15 


0.0089 


0.0009 


9.05 


0.31 


0.3238 


1 


H 


RXJ1756.1+6615 


160 


17 56 10.0 


+66 15 


14 


17 56 09.5 


+66 


15 


09 


4.01 


0.0086 


0.0007 


8.60 


1.48 


0.6357 


1 


H 
d 


RX .11756.1+7001 


3460 


17 56 10.8 


+70 01 


55 


17 56 11.9 


+70 


01 


47 


3.90 


0.0090 


0.0017 


8.88 


0.12 


0.2129 


1 


RX .11756.1+7055 


3470 


17 56 10.8 


+70 55 


48 


17 56 11.6 


+70 


55 


50 


4.25 


0.0163 


0.0022 


16.76 


0.31 


0.2460 


1 


> 


RXJ1756.2+6619 


170 


17 56 12.0 


+66 19 


47 


17 56 12.2 


+66 


19 


46 


4.01 


0.0023 


0.0004 


2.30 


1.64 


1.1340 


1 


O 


RXJ1756.2+6955 


3480 


17 56 12.7 


+69 55 


21 


17 56 12.7 


+69 


55 


20 


3.84 


0.0177 


0.0023 


17.33 


0.03 


0.0838 


2 


RX J1756.2+6838 


3490 


17 56 13.6 


+68 38 


31 


17 56 15.5 


+68 


38 


25 


5.24 


0.0070 


0.0012 


7.92 


0.02 


0.1019 


2 




RX J1756.4+6300 


3510 


17 56 25.2 


+63 00 


42 


17 56 25.4 


+63 


00 


49 


3.53 


0.0077 


0.0016 


7.23 


4.90 


1.1110 


1 




RXJ1756.7+6438 


3530 


17 56 43.2 


+64 38 


53 


17 56 43.4 


+64 


38 


59 


3.49 


0.0231 


0.0019 


21.57 


0.32 


0.2233 


1 




RX J1756.8+6612 


180 


17 56 52.4 


+66 12 


42 


17 56 51.3 


+66 


12 


42 


4.14 


0.0061 


0.0006 


6.20 


7.78 


1.4252 


1 




RX J1756.9+6238 


3550 


17 56 58.2 


+62 38 


44 


17 57 00.6 


+62 


38 


55 


3.33 


0.0205 


0.0026 


18.70 


12.10 


1.0902 


1 




RXJ1757. 1+6352 


3570 


17 57 09.7 


+63 52 


38 


17 57 09.5 


+63 


52 


33 


3.31 


0.0145 


0.0019 


13.19 


0.45 


0.3220 


1 




RXJ1757.5+6841 


3590 


17 57 34.1 


+68 41 


22 


17 57 34.1 


1 go 
-t-0c5 




21 


5.46 


0.0205 


0.0018 


23.61 


0.22 


0.1814 


1 




RX .71757.9+6934 


3600 


17 57 55.2 


+69 34 


23 


17 57 55.2 


+69 


34 


25 


4.35 


0.0270 


0.0024 


28.09 


0.04 


0.0795 


1 




RX.) 1757.9+6609 


210 


17 57 56.9 


+66 09 


23 


17 57 56.5 


+66 


09 


20 


4.26 


0.0025 


0.0005 


2.58 


0.23 


0.4865 


2 




RXJ1758.0+6851 


3620 


17 58 02.4 


+68 51 


46 


17 58 03.7 


+68 


51 


51 


4.87 


0.0061 


0.0013 


6.69 


0.07 


0.1876 


1 




RXJ1758.2+7020 


3621 


17 58 12.4 


+70 20 


27 


17 58 13.4 


+70 


20 


23 


4.25 


0.0080 


0.0018 


8.23 


5.94 


1.1400 


1 




RXJ1758.2+6743 


230 


17 58 13.2 


+67 43 


18 


17 58 14.1 


+67 


43 


17 


5.02 


0.0122 


0.0012 


13.55 


0.16 


0.2045 


1 




RXJ1758.3+6906 


3630 


17 58 18.7 


+69 06 


30 


17 58 15.9 


+69 


06 


32 


4.42 


0.0104 


0.0016 


10.91 


37.92 


2.1572 


1 




RXJ1758.3+6203 


241 


17 58 23.3 


+62 03 


26 


17 58 24.4 


+62 


03 


20 


3.41 


0.0078 


0.0019 


7.21 


1.35 


0.6590 


1 




RX .11758.4+6531 


250 


17 58 24.1 


+65 31 


05 


17 58 24.2 


+65 


31 


08 


4.04 


0.0195 


0.0013 


19.57 


0.68 


0.3250 


1 




RXJ1758.7+6423 


3680 


17 58 44.5 


+64 23 


04 


17 58 43.1 


+64 


23 


04 


3.38 


0.0054 


0.0011 


4.97 


1.29 


0.7523 


1 




RXJ1758.8+6551 


280 


17 58 52.8 


+65 51 


06 


17 58 53.2 


+65 


51 


13 


4.40 


0.0030 


0.0005 


3.14 


0.17 


0.3884 


2 




RXJ1758.9+6220 


3700 


17 58 56.5 


+62 20 


31 


17 58 56.3 


+62 


20 


25 


3.33 


0.0070 


0.0017 


6.38 


1.34 


0.6910 


1 


h- ' 
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03 



Object 




NEP 




ax 


<5x 




aopt 


5, 








Count Rate 


Count Rate 


/x 


ix 


z 


AGN 








ID 


(J2000) 


(J2000) 


(J2000) 


(J2000) 


[10^0 cm-2] 


[s-l] 


Error [s" 


[10-1" cgs] 


[1044 cgs] 




Type 




(1) 




(2) 




(3) 


(4) 




(5) 


(6) 




(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 




RXJ1759.34 


-6335 


3730 


17 


59 


19.2 


+63 35 37 


17 


59 


18.4 


+63 


35 


40 


3.15 


0.0063 


0.0014 


5.59 


4.92 


1.2354 


1 




RXJ1759.7+6739 


350 


17 


59 


42.5 


+67 39 25 


17 


59 


42.8 


+67 


39 


27 


5.35 


0.0049 


0.0007 


5.60 


3.56 


1.0830 


1 




RXJ1759.7+6629 


360 


17 


59 


44.3 


+66 29 11 


17 


59 


44.7 


+66 


29 


11 


4.12 


0.0040 


0.0004 


4.06 


0.23 


0.3990 


1 




RXJ1800. 1+6636 


380 


18 


00 


07.5 


+66 36 54 


18 


00 


07.6 


+66 


36 


55 


4.17 


0.0048 


0.0004 


4.89 


0.001 


0.0260 


2 




RXJ1800. 1+6938 


381 


18 


00 


08.4 


+69 38 30 


18 


00 


10.4 


+69 


38 


38 


4.56 


0.0055 


0.0013 


5.85 


3.57 


1.0650 


1 




RXJ1800. 1+6720 


390 


18 


00 


11.2 


+67 20 48 


18 


00 


10.9 


+67 


20 


58 


4.68 


0.0036 


0.0006 


3.87 


2.82 


1.1433 


1 




RXJ1800.3+6615 


430 


18 


00 


23.1 


+66 15 54 


18 


00 


23.8 


+66 


15 


52 


4.03 


0.0048 


0.0005 


4.81 


0.35 


0.4475 


1 




RXJ1800.4+7051 


3810 


18 


00 


25.2 


+70 51 55 


18 


00 


25.8 


+70 


51 


58 


4.25 


0.0101 


0.0019 


10.39 


0.35 


0.3200 


1 




RXJ1800.4+6357 


3820 


18 


00 


26.2 


+63 57 19 


18 


00 


26.4 


+63 


57 


19 


2.72 


0.0205 


0.0020 


16.83 


3.44 


0.6828 


1 




RXJ1800.44 


-6705 


440 


18 


00 


29.0 


+67 05 48 


18 


00 


28.9 


+67 


05 


50 


4.52 


0.0041 


0.0005 


4.34 


4.62 


1.3330 


1 




RXJ1801.24 


-6433 


3860 


18 


01 


13.2 


+64 33 22 


18 


01 


14.7 


+64 


33 


26 


3.43 


0.0061 


0.0011 


5.65 


2.09 


0.8700 


1 




RXJ1801.24 


-6902 


3870 


18 


01 


14.6 


+69 02 43 


18 


01 


14.6 


+69 


02 


43 


4.86 


0.0057 


0.0012 


6.24 


5.89 


1.2700 


1 




RXJ1801.24 


-6624 


480 


18 


01 


15.2 


+66 24 01 


18 


01 


16.6 


+66 


24 


01 


4.05 


0.0012 


0.0002 


1.21 


1.09 


1.2500 


1 




RXJ1802.04 


-6629 


560 


18 


02 


05.9 


+66 29 02 


18 


02 


04.8 


+66 


29 


14 


4.14 


0.0013 


0.0003 


1.32 


0.03 


0.2650 


1 




RXJ1802.14 


-6535 


570 


18 


02 


07.7 


+65 35 21 


18 


02 


07.7 


+65 


35 


14 


4.23 


0.0034 


0.0007 


3.49 


0.02 


0.1513 


1 




RXJ1802.34 


-6259 


3910 


18 


02 


19.6 


+62 59 21 


18 


02 


21.5 


+62 


59 


14 


3.48 


0.0063 


0.0014 


5.88 


1.39 


0.7240 


1 




RXJ1802.34 


-6647 


590 


18 


02 


22.8 


+66 47 49 


18 


02 


24.5 


+66 


47 


35 


4.24 


0.0047 


0.0005 


4.83 


0.19 


0.3424 


1 




RXJ1802.74 


-6727 


630 


18 


02 


47.4 


+67 27 50 


18 


02 


47.8 


+67 


27 


41 


4.94 


0.0026 


0.0006 


2.87 


0.02 


0.1620 


2 




RX .11802.84 


-6605 


640 


18 


02 


51.2 


+66 05 40 


18 


02 


51.3 


+66 


05 


42 


4.12 


0.0166 


r\ c\r\c\c\ 

0.U009 


16.82 


0.21 


0.2070 


1 




RX J 1803.4+6738 


650 


18 


03 


28.3 


+67 38 06 


18 


03 


29.0 


+67 


38 


10 


4.87 


0.2801 


0.0042 


306.77 


1.51 


0.1360 


1 




RXJ1803.8+6619 


660 


18 


03 


50.4 


+66 19 31 


18 


03 


50.1 


+66 


19 


31 


4.30 


0.0064 


0.0006 


6.62 


0.37 


0.3968 


1 




RXJ1804.5+6937 


4020 


18 


04 


34.2 


+69 37 33 


18 


04 


34.4 


+69 


37 


37 


4.76 


0.0082 


0.0016 


8.90 


1.36 


0.6055 


1 


ci 


RXJ1804.6+6846 


4040 


18 


04 


41.9 


+68 46 02 


18 


04 


40.6 


+68 


45 


55 


5.07 


0.0072 


0.0013 


8.03 


0.02 


0.0969 


2 




RXJ1805.2+7006 


4080 


18 


05 


16.6 


+70 06 19 


18 


05 


17.8 


+70 


06 


22 


4.19 


0.0162 


0.0024 


16.55 


0.16 


0.1874 


1 




RXJ1805.4+6638 


710 


18 


05 


25.3 


+66 38 58 


18 


05 


25.0 


+66 


39 


03 


4.10 


0.0105 


0.0008 


10.61 


0.06 


0.1449 


1 


1— 1 

m 


RXJ1805.6+6624 


720 


18 


05 


36.1 


+66 24 52 


18 


05 


36.2 


+66 


24 


52 


4.10 


0.0110 


0.0008 


11.12 


2.60 


0.7210 


1 


H 


RXJ1805.6+6309 


4110 


18 


05 


39.0 


+63 09 36 


18 


05 


40.1 


+63 


09 


22 


3.00 


0.0080 


0.0016 


6.92 


0.67 


0.5013 


1 




RXJ1805.64 


-6432 


4120 


18 


05 


41.4 


+64 32 51 


18 


05 


40.5 


+64 


32 


47 


3.35 


0.0121 


0.0015 


11.07 


2.78 


0.7432 


1 


> 


RXJ1806.04 


-6940 


4140 


18 


06 


03.2 


+69 40 26 


18 


06 


03.2 


+69 


40 


24 


4.76 


0.0087 


0.0017 


9.44 


0.32 


0.3214 


1 




RXJ1806.24 


-6644 


740 


18 


06 


12.5 


+66 44 40 


18 


06 


12.4 


+66 


44 


34 


4.00 


0.0041 


0.0006 


4.10 


0.17 


0.3482 


1 




RXJ1807.74 


-6617 


4270 


18 


07 


47.4 


+66 17 32 


18 


07 


47.4 


+66 


17 


31 


3.94 


0.0087 


0.0008 


8.63 


3.82 


0.9350 


1 




RX .11808.04 


-6452 


4280 


18 


08 


02.5 


+64 52 24 


18 


08 


03.7 


+64 


52 


30 


3.38 


0.0211 


0.0021 


19.38 


11.05 


1.0360 


1 




RXJ1808.84 


-6634 


4390 


18 


08 


49.8 


+66 34 31 


18 


08 


49.6 


+66 


34 


29 


3.85 


0.0220 


0.0013 


21.59 


4.64 


0.6970 


1 




RXJ1808.84 


-6530 


4400 


18 


08 


51.0 


+65 30 21 


18 


08 


50.8 


+65 


30 


19 


3.91 


0.0095 


0.0013 


9.38 


0.26 


0.2937 


2 




RXJ1808.84 


-6511 


4401 


18 


08 


53.4 


+65 11 42 


18 


08 


53.5 


+65 


11 


48 


3.82 


0.0043 


0.0010 


4.20 


7.40 


1.6350 


1 




RXJ1809.04 


-6704 


4410 


18 


09 


01.0 


+67 04 21 


18 


09 


00.9 


+67 


04 


25 


4.09 


0.0046 


0.0008 


4.65 


0.99 


0.6950 


1 




RXJ1809.0+6800 


4420 


18 


09 


03.5 


+68 00 55 


18 


09 


03.1 


+68 


00 


57 


4.47 


0.0062 


0.0012 


6.53 


0.95 


0.5946 


1 




RXJ1809.0+6333 


4430 


18 


09 


05.0 


+63 33 00 


18 


09 


05.5 


+63 


33 


09 


2.65 


0.0073 


0.0017 


5.92 


1.04 


0.6412 


1 




RXJ1809.5+6620 


4440 


18 


09 


30.1 


+66 20 33 


18 


09 


30.1 


+66 


20 


21 


3.86 


0.0061 


0.0009 


5.99 


1.03 


0.6350 


1 




RXJ1809.5+6609 


4450 


18 


09 


34.8 


+66 09 06 


18 


09 


34.2 


+66 


09 


11 


3.85 


0.0055 


0.0009 


5.40 


2.42 


0.9400 


1 




RXJ1809.7+6837 


4460 


18 


09 


46.8 


+68 37 26 


18 


09 


48.4 


+68 


37 


34 


5.36 


0.0059 


0.0012 


6.74 


0.09 


0.2173 


1 




RXJ1810.0+6344 


4490 


18 


10 


04.2 


+63 44 24 


18 


10 


04.4 


+63 


44 


26 


2.76 


0.0223 


0.0026 


18.47 


0.91 


0.3770 


1 




RXJ1810.3+6328 


4501 


18 




Zo.O 


+63 28 08 


18 


1 fi 


io.y 


+0o 




14 


2.74 


0.0135 


0.0025 


11.12 


3.76 


0.8380 


1 




RXJ1810.4+6432 


4520 


18 


10 


24.7 


+64 32 46 


18 


10 


24.2 


+64 


32 


54 


3.12 


0.0067 


0.0015 


5.91 


0.17 


0.3030 


1 




RXJ1811. 2+6543 


4550 


18 


11 


12.4 


+65 43 46 


18 


11 


11.6 


+65 


43 


47 


3.94 


0.0127 


0.0014 


12.59 


1.15 


0.4895 


1 




RXJ1811.64 


-6507 


4580 


18 


11 


36.8 


+65 07 04 


18 


11 


36.1 


+65 


07 


00 


3.67 


0.0197 


0.0020 


18.86 


6.54 


0.8470 


1 




RXJ1811.64 


-6333 


4oyu 


18 


11 


41.2 


1 fiq QQ AP. 


18 


11 


43.5 


+63 


33 


51 


2 72 




n nni Q 
u.uuiy 


u.yu 




U.OOIU 


1 




RXJ1812.44 


-6610 


4640 


18 


12 


27.0 


+66 10 46 


18 


12 


26.1 


+66 


10 


48 


3.87 


0.0036 


0.0008 


3.54 


0.63 


0.6449 


1 




RXJ1813.04 


-6644 


4670 


18 


13 


04.8 


+66 44 56 


18 


13 


06.1 


+66 


44 


52 


4.56 


0.0056 


0.0010 


5.96 


1.62 


0.7680 


1 




RXJ1813.14 


-6547 


4680 


18 


13 


09.0 


+65 47 01 


18 


13 


07.7 


+65 


47 


04 


3.87 


0.0102 


0.0013 


10.04 


0.41 


0.3489 


1 




RXJ1813.14 


-6608 


4690 


18 


13 


10.7 


+66 08 02 


18 


13 


07.9 


+66 


08 


09 


3.83 


0.0045 


0.0009 


4.40 


4.75 


1.3400 


1 




RXJ1813.54 


-6635 


4720 


18 


13 


34.1 


+66 35 36 


18 


13 


35.1 


+66 


35 


34 


4.58 


0.0038 


0.0008 


4.05 


0.76 


0.6609 


1 




RXJ1813.64 


-6731 


4721 


18 


13 


41.5 


+67 31 50 


18 


13 


43.0 


+67 


32 


23 


4.67 


0.0142 


0.0019 


15.27 


2.44 


0.6168 


1 




RXJ1813.7+6538 


4760 


18 


13 


46.6 


+65 38 21 


18 


13 


45.8 


+65 


38 


20 


3.91 


0.0342 


0.0022 


33.79 


0.35 


0.1912 


1 
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Object NEP ax <5x aopt 5opt nn Count Rate Count Rate /x ^x ^ AGN 
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(1) 




\'^) 












{0) 


(6) 




{' ) 


^^^^ 






n 1 "i 


(LZ) 






RXJ1813.8+6728 


4800 


18 


13 51.0 


+67 28 


10 


18 


13 50.7 


+67 


28 


06 


4.63 


0.0084 


0.0015 


9.00 


0.30 


0.3196 


1 




RXJ1815.2+6658 


4880 


18 


15 17.0 


+66 58 


10 


18 


15 17.4 


+66 


58 


05 


4.67 


0.0070 


0.0012 


7.53 


0.12 


0.2287 


2 




RXJ1815.3+6507 


4890 


18 


15 19.1 


+65 07 


28 


18 


15 20.0 


+65 


07 


14 


3.85 


0.0073 


0.0015 


7.16 


14.36 


1.7234 


1 




RXJ1815.4+6806 


4910 


18 


15 24.4 


+68 06 


29 


18 


15 24.9 


+68 


06 


32 


4.58 


0.0214 


0.0022 


22.82 


0.39 


0.2390 


1 




RX J1815.8+6441 


4930 


18 


15 52.3 


+64 41 


00 


18 


15 51.7 


+64 


41 


03 


3.16 


0.0117 


0.0020 


10.38 


0.63 


0.4116 


1 




RXJ1817.5+6631 


5020 


18 


17 32.1 


+66 31 


08 


18 


17 31.3 


+66 


31 


11 


4.64 


0.0083 


0.0012 


8.90 


0.67 


0.4513 


1 




RXJ1818.4+6741 


5050 


18 


18 28.9 


+67 41 


26 


18 


18 28.8 


+67 


41 


24 


4.64 


0.1125 


0.0044 


120.64 


3.87 


0.3140 


1 




RXJ1818.7+6518 


5080 


18 


18 46.2 


+65 18 


14 


18 


18 44.8 


+65 


18 


10 


4.20 


0.0056 


0.0013 


5.73 


5.82 


1.3080 


1 




RXJ1819.8+6510 


5150 


18 


19 52.2 


+65 10 


35 


18 


19 51.5 


+65 


10 


37 


4.21 


0.0173 


0.0020 


17.72 


0.18 


0.1894 


1 




RXJ1819.9+6628 


5190 


18 


19 59.9 


+66 28 


25 


18 


19 59.0 


+66 


28 


30 


4.79 


0.0049 


0.0012 


5.33 


5.07 


1.2740 


1 




RXJ1820.5- 


h6620 


5230 


18 


20 32.9 


+66 20 


29 


18 


20 ,32.9 


+66 


20 


20 


4.71 


0.0075 


0.0015 


8.09 


0.80 


0.5057 


1 




RXJ1820.5- 


h6930 


5240 


18 


20 35.3 


+69 30 


04 


18 


20 ,33.4 


+69 


30 


15 


6.46 


0.0160 


0.0025 


19.77 


0.24 


0.2051 


1 


o 


RXJ1821.6- 


h6543 


5300 


18 


21 38.8 


+65 43 


04 


18 


21 40.0 


+65 


43 


10 


4.26 


0.0171 


0.0020 


17.63 


0.39 


0.2666 


2 


RXJ1821.6- 


h6328 


5310 


18 


21 39.6 


+63 28 


27 


18 


21 38.8 


+63 


28 


26 


3.10 


0.0282 


0.0035 


24.80 


1.13 


0.3656 


1 


a 


RXJ1821.9- 


h6654 


5321 


18 


21 55.7 


+66 54 


34 


18 


21 56.5 


+66 


54 


26 


5.07 


0.0056 


0.0011 


6.24 


0.01 


0.0873 


1 




RXJ1821.9- 


h6420 


5340 


18 


21 57.4 


+64 20 


51 


18 


21 57.1 


+64 


20 


37 


3.77 


1.0710 


0.0140 


1039.63 


29.32 


0.2970 


1 


lERING 


RXJ1821.9- 


h6818 


5330 


18 


21 58.8 


+68 18 


42 


18 


21 59.4 


+68 


18 


42 


4.85 


0.0073 


0.0017 


7.98 


15.30 


1.6920 


1 


RXJ1823.3" 


h6419 


5400 


18 


23 20.0 


+64 19 


23 


18 


23 19.2 


+64 


19 


32 


3.32 


0.0219 


0.0026 


19.95 


2.70 


0.5766 


1 


RXJ1823.6- 


h6847 


5411 


18 


23 38.6 


+68 47 


40 


18 


23 39.4 


+68 


47 


46 


5.85 


0.0085 


0.0019 


10.08 


0.13 


0.2071 


1 


RXJ1823.9+6719 


5440 


18 


23 54.6 


+67 19 


41 


18 


23 54.7 


+67 


19 


36 


4.75 


0.0079 


0.0015 


8.56 


0.65 


0.4536 


1 


RXJ1824.7+6509* 


5500 


18 


24 46.9 


+65 09 


24 


18 


24 47.3 


+65 


09 


25 


3.97 


0.1269 


0.0048 


126.39 


3.73 


0.3030 


1 


O 


RXJ1825.7+6905 


5530 


18 


25 46.3 


+69 05 


51 


18 


25 47.3 


+69 


05 


54 


6.66 


0.0362 


0.0036 


45.27 


0.09 


0.0888 


1 




RXJ1826.6+6706 


5550 


18 


26 38.3 


+67 06 


47 


18 


26 37.5 


+67 


06 


44 


5.48 


0.0212 


0.0023 


24.46 


0.64 


0.2870 


1 




RXJ1827.2+6549 


5560 


18 


27 15.3 


+65 49 


21 


18 


27 13.9 


+65 


49 


20 


5.30 


0.0071 


0.0016 


8.07 


6.97 


1.2250 


1 


pi 


RXJ1827.5+6431 


5590 


18 


27 33.6 


+64 31 


38 


18 


27 33.8 


+64 


31 


44 


4.16 


0.0141 


0.0022 


14.35 


0.03 


0.0977 


1 




RXJ1828. 1+6709 


5601 


18 


28 06.6 


+67 09 


23 


18 


28 06.7 


+67 


09 


17 


6.02 


0.0085 


0.0016 


10.20 


4.61 


0.9430 


1 




RXJ1828.2+6403 


5620 


18 


28 13.7 


+64 03 


31 


18 


28 14.2 


+64 


03 


28 


3.83 


0.0119 


0.0023 


11.64 


0.03 


0.0963 


1 




RXJ1828.7- 


h6953 


5670 


18 


28 47.9 


+69 53 


58 


18 


28 49.3 


+69 


54 


00 


7.02 


0.0271 


0.0033 


34.60 


0.11 


0.1100 


1 


H 


RXJ1828.8- 


h6452 


5680 


18 


28 48.6 


+64 52 


50 


18 


28 48.3 


+64 


53 


00 


4.42 


0.0134 


0.0020 


14.05 


5.25 


0.8730 


1 


t-i 
H 


RXJ1829.0- 


h6433 


5690 


18 


29 00.6 


+64 33 


49 


18 


29 00.4 


+64 


33 


51 


4.30 


0.0183 


0.0023 


18.94 


1.00 


0.3880 


1 


O 


RXJ1829.5- 


h6631 


5740 


18 


29 35.4 


+66 31 


19 


18 


29 35.4 


+66 


31 


23 


6.21 


0.0066 


0.0016 


8.03 


1.68 


0.6898 


1 


H 


RXJ1829.7- 


h6749 


5750 


18 


29 43.4 


+67 49 


09 


18 


29 42.2 


+67 


49 


12 


6.19 


0.0158 


0.0023 


19.18 


1.66 


0.4783 


1 


H 
d 


RXJ1830.0- 


h6645 


5790 


18 


30 01.4 


+66 45 


23 


18 


,30 02.0 


+66 


45 


23 


6.38 


0.0142 


0.0021 


17.45 


0.46 


0.2889 


1 


RXJ1830.1- 


h642.5 


5800 


18 


30 07.5 


+64 25 


28 


18 


,30 06.1 


+64 


25 


29 


3.90 


0.0104 


0.0019 


10.26 


1.69 


0.6253 


1 


> 


RXJ1832.0- 


h6542 


5890 


18 


32 01.5 


+65 42 


35 


18 


32 01.3 


+65 


42 


33 


5.41 


0.0239 


0.0026 


27.42 


0.74 


0.2908 


1 


O 


RXJ1832.0- 


h6447 


5900 


18 


32 04.2 


+64 47 


01 


18 


32 01.2 


+64 


47 


08 


4.33 


0.0177 


0.0025 


18.37 


12.65 


1.1180 


1 


RXJ1832.4+6402 


5930 


18 


32 25.2 


+64 02 


02 


18 


32 24.1 


+64 


02 


10 


4.93 


0.0106 


0.0022 


11.68 


0.59 


0.3826 


1 




RXJ1832.4+6438 


5940 


18 


32 25.2 


+64 38 


15 


18 


32 24.2 


+64 


38 


23 


4.27 


0.0140 


0.0021 


14.45 


1.62 


0.5335 


1 




RXJ1833.0+6344 


5990 


18 


33 02.8 


+63 44 


17 


18 


33 01.3 


+63 


44 


35 


4.93 


0.0139 


0.0024 


15.31 


1.88 


0.5535 


1 




RXJ1835.0+6526 


6070 


18 


35 04.8 


+65 26 


44 


18 


35 06.0 


+65 


26 


47 


5.52 


0.0206 


0.0026 


23.85 


1.41 


0.4083 


1 




RXJ1835. 1+6342 


6080 


18 


35 08.0 


+63 42 


33 


18 


35 10.0 


+63 


43 


14 


4.93 


0.0228 


0.0033 


25.12 


11.40 


0.9445 


1 




RXJ1835. 1+6733 


6090 


18 


35 10.3 


+67 33 


54 


18 


35 09.0 


+67 


33 


58 


6.59 


0.0075 


0.0017 


9.34 


2.37 


0.7460 


1 




RXJ1836.4+6602 


6180 


18 


36 28.2 


+66 02 


40 


18 


36 28.7 


-f-OD 


no 


37 


7.08 


0.0130 


0.0024 


16.65 


0.16 


0.1858 


1 




RX .11836.6+6719 


6200 


18 


36 36.1 


+67 19 


04 


18 


36 ,36.8 


+67 


19 


12 


8.03 


0.0084 


0.0018 


11.29 


0.25 


0.2693 


1 




RX .11838.1+6649 


6280 


18 


38 09.1 


+66 49 


26 


18 


38 10.0 


+66 


49 


22 


7.17 


0.0083 


0.0019 


10.69 


0.28 


0.2879 


1 




RXJ1838.8- 


h6432 


6300 


18 


38 51.9 


+64 32 


21 


18 


,38 53.1 


+64 


32 


23 


5.43 


0.0081 


0.0019 


9.31 


0.77 


0.4700 


1 




RXJ1839.2- 


h6711 


6301 


18 


39 16.9 


+67 11 


12 


18 


,39 16.5 


+67 


11 


06 


8.24 


0.0063 


0.0015 


8.55 


1.94 


0.7130 


1 




RXJ1839.3- 


h6544 


6340 


18 


39 18.5 


+65 44 


42 


18 


,39 18.3 


+65 


44 


35 


5.77 


0.0095 


0.0023 


11.20 


0.02 


0.0820 


1 




RX .11841. 3- 


h6321 


6450 


18 


41 18.9 


+63 21 


36 


18 


41 20.0 


+63 


21 


42 


5.32 


0.0096 


0.0022 


10.94 


15.55 


1.4990 


1 




RXJ1842.2- 


h6204 


6452 


18 


42 14.8 


+62 04 


24 


18 


42 15.6 


+62 


04 


24 


5.25 


0.0146 


0.0032 


16.54 


0.56 


0.3203 


1 




RXJ1842.5- 


h6809 


6490 


18 


42 33.0 


+68 09 


30 


18 


42 33.3 


+68 


09 


25 


6.35 


0.0511 


0.0044 


62.70 


5.33 


0.4750 


1 




RXJ1842.9- 


h6241 


6491 


18 


42 56.4 


+62 41 


44 


18 


42 55.2 


+62 


41 


49 


5.32 


0.0502 


0.0048 


57.21 


0.10 


0.0835 


1 




RXJ1843.3- 


h6653 


6520 


18 


43 22.5 


+66 53 


21 


18 


43 20.9 


+66 


53 


29 


7.33 


0.0110 


0.0022 


14.28 


0.50 


0.3273 


1 




RXJ1843.9+6821 


6540 


18 


43 55.7 


+68 21 


11 


18 


43 54.1 


+68 


21 


01 


6.08 


0.0200 


0.0032 


24.10 


1.12 


0.3688 


1 


h- ' 



TABLE 1 — Continued 



00 



Object 


NEP 


ax 


<5x 


Qopt 


<5opt 


riH 


Count Rate 


Count Rate 


/x 




z 


AGN 




ID 


(J2000) 


(J2000) 


(J2000) 


(J2000) 


[lO^" cm-2] 


[s-i] 


Error 


[10"" cgs] 


[1044 cgs] 




Type 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


RXJ1844.3+6431 


6542 


18 44 23.4 


+64 31 31 


18 44 21.8 


+64 31 46 


6.00 


0.0105 


0.0022 


12.59 


1.13 


0.4870 


1 


RXJ1844.4+6236 


6543 


18 44 26.9 


+62 36 12 


18 44 26.4 


+62 36 14 


5.32 


0.0274 


0.0039 


31.22 


1.03 


0.3172 


1 


RXJ1844.4+6248 


6544 


18 44 27.5 


+62 48 27 


18 44 26.2 


+62 48 29 


5.32 


0.0239 


0.0034 


27.24 


67.62 


1.8800 


1 


RXJ1844.9+6813 


6570 


18 44 54.0 


+68 13 23 


18 44 54.1 


+68 13 17 


5.84 


0.0250 


0.0036 


29.64 


0.92 


0.3097 


1 



* RX J1824. 7+6509 (NEP 5500) was formerly classified as a stellar X-ray source based on the literature. However, more recent work convincingly demonstrates this source is a type 1 AGN at 
z = 0.303 fe.g-. IEngels et aHl 998). 

Note — Table 1 is also available in machine-readable form in the eletronic edition of the Astrophysical Journal. 
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